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Abstract 


In this paper, the theory of plausible and paradoxical reasoning of Dezert- 
Smarandache (DSmT) is used to take into account the paradoxical character 
through the intersections of vegetation, aquatic and mineral surfaces. In order 
to do this, we developed a classification model of pixels by aggregating informa- 
tion using the DSmT theory based on the PCR5 rule using the NDVI, MNDW1 
and NDBal spectral indices obtained from the ASTER satellite images. On the 
qualitative level, the model produced three simple classes for certain knowledge 
(E, V, M) and eight composite classes including two union classes characterizing 
partial ignorance ({E,V}, {M,V}) and six classes of intersection of which three 
classes of simple intersection (ENV, MAV, EOM) and three classes of compo- 
site intersection (EX{M,V}, MO{E,V}, VA{E,M}), which represent paradoxes. This 
model was validated with an average rate of 93.34% for the well-classified pix- 
els and a compliance rate of the entities in the field of 96.37%. Thus, the model 1 
retained provides 84.98% for the simple classes against 15.02% for the composite 
classes. 


Keywords 
Theory of Plausible and Paradoxical Reasoning, PCR5 Rule, ASTER Satellite 
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Images, NDVI, MNDWI, NDBal, Modeling, Classification 


1. Introduction 


The surface state of the Earth can be represented in remote sensing by three ent- 
ities that are vegetation surfaces, aquatic surfaces and mineral surfaces. Any sur- 
face area observed, depending on the size of the area, may be a combination of 
these three entities. 

Thus, it can be observed vegetation surface, aquatic surface, mineral surface, 
vegetation and aquatic surface, vegetation and mineral surface, mineral and aq- 
uatic surface, vegetation, mineral and aquatic surface. 

The use of satellite images, for mapping purposes, has been the subject of sev- 
eral studies ([1] [2] [3] [4]). Researchers have used spectral indices to map vege- 
tation surfaces, aquatic surfaces, and surfaces of bare soil and mansions. 

These indices were used for image classification. Unfortunately, imperfections 
(uncertainties, inaccuracies, etc.) on the information produced by the images 
associated with said indices are observed ([5] [6]). Taking into account and good 
management of these imperfections are done by the theory of belief functions. 

The Dempster-Shafer theory has been used successfully to handle cases of un- 
certainty, vagueness and ignorance in the classification of pixels to classes of ve- 
getation surface, aquatic surface and mineral surface [4]. 

In this article, the theory of plausible and paradoxical reasoning of Dezert- 
Smarandache (DSmT) is used to take into account the paradoxical character through 
the intersections of the elements vegetation surface, aquatic surface and mineral 
surface. 

The general objective of the study is to develop a model of pixel classification 
by aggregating information using the DSmT theory, spectral indices NDVI (Nor- 
malized Difference Vegetation Index), MNDWI (Modification of Normalized Dif- 
ference Water Index) and NDBa/ (Normalized Difference Bare Index) and ASTER 
satellite images. It acts specifically first, to model the frameworks of discernment 
and reasoning and belief functions, then define the decision criteria and write 
algorithms and programming codes under the MATLAB software; finally realize 
and evaluate classified image. 

The contribution of this study is to give an approach of unsupervised classifi- 
cation of mapping that takes into account the plausible and paradoxical charac- 
teristics related to the information of ground that is to say without needing to 
know the real spatial state of the ground concerned. 

This paper, which proposes to report on the work carried out, presents succes- 
sively the theory of plausible and paradoxical reasoning of Dézert-Smarandache, 
the material used, the methodological approach that guided the work and the re- 


sults obtained. 





DOI: 10.4236/ojapps.2017.710038 


521 Open Journal of Applied Sciences 


J.-C. Okaingni et al. 


2. Principle and Formalism of DSmT 


Dezert-Smarandache theory (DSmT) is interpreted as a generalization of the 
theory of Dempster-Shafer (DST) ([7]). Its basic principle and formalism for ag- 
gregating information can be characterized by a four-stage structure of modeling, 
estimation, combination and decision. 

Modeling consists in choosing the representation of the frameworks of dis- 
cernment and reasoning and the models of the mass functions to be used. 

The discernment framework © ={6,,---,0,} of DSmT is an exhaustive set of 
different hypotheses, not necessarily exclusive assumptions. The reasoning frame- 
work D? associated, also called hyper-powerset, is the set of all possible proposi- 
tions constructed from the elements of ©, including the empty set (ø), with the 
operators U and A. It is characterized by the following conditions ([8]): 

1) 9,6,,°:-,Oy € De; 

2)If A,A €D° then A,X A,eD®° and A,UA,€D*, 

3) There are no other elements belonging to D®, except those obtained using rules 
1 or 2. 

The construction of the reasoning framework D? can be obtained by a matrix 
product between the binary matrix of Dedekind and the coding vector of Sma- 
randache ([1] [7]). Its cardinal increases according to the cardinal of the discern- 
ment framework on which it is based. 

In order to fix ideas, it is considered the cardinal of © equal to 3. Thus, ac- 
cording to Dedekind ((9]), the cardinal of D® is equal to 19. Then we get ([1]): 


© ={6,,0,,0,} 
D° = {4,0,,0,,0;,0, J ,0, 90,0, VA, 0, VO, V8,,8,.08;, 
9, 09;,8,08;,8, 1(8, VO; ),A; (8, VA, ),8, A(8, U4), 
8, (8, 08;),O; U(8,09,),4, (8.08), 
(6,18,)U(8,. 18, )U(8, 08),8, 09, 08,} 


In general, we define a generalized mass function m of DP with values in [0,1] 


satisfying the following conditions of Equation (1): 


£ m(A)=1 
AeD® (1) 
m(¢@) =0 


where ¢ is the empty set. The value mA) quantifies the belief that the class 
sought belongs to the subset A of D? (and to no other subset of A). The subsets 
A such that m(A) > 0 are called focal elements. 
The following special mass functions are defined ([10]): 
v amass function mis said to be normal when m(¢) = 0; 
v amass function m is said to be categorical when it has a single focal element 
A such that m(A) = 1. In the case where A is a set, knowledge is certain but 
imprecise. 


When A = {G}, knowledge is certain and precise; 
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v amass function m is said to be empty (or total ignorance) when the func- 
tion m is categorical in ©: m(@) = 1; 
v amass function m is said to be Bayesian if all the focal elements are single- 
tons of © in Equation (2): 
> m(0,)=1; (2) 
Q; €O 
v amass function m is said to be dogmatic if m(®) = 0; 
v amass function zm is said to be consonant if all the focal elements are nested; 
v amass function m is said to have a simple support when it has 2 focal ele- 


ments, one of which is © (Equation (3)): 
m(A) =l-w, AeD? 
m(©) =O (3) 
m(B)=0, VBeD*°\A 


In this case, the function m can also be denoted A® where w represents the 
weight of the ignorance of the mass function m. 

The estimation consists in determining all the parameters of the mass func- 
tions selected at the modeling stage. This is a difficult problem that does not have a 
universal solution. The difficulty is further increased if we want to assign masses 
to compound hypotheses involving intersections and/or unions ([11]). 

The combination is the grouping phase of the information, from the mass func- 
tions of the different information sources, using an operator adapted to the for- 
malism of the modeling. The DSmT has two types of combination ([12]): the clas- 
sic version and the hybrid version. Consider n initial mass functions m,,m,,---,m, 
representing the respective information of n different sources, which can be com- 
bined according to the DSmT. 

The classic combination of DSmT is (Equation (4)): 

m(C)= ` m (A, )*m, (A, )*---*m, (A, ) 
AYO AQ OA, =C (4) 


A,(i=1---,n)€D°, VC € D° 


The hybrid combination of DSmT is used in the presence of integrity con- 
straints applied to D®. An integrity constraint of a set U is an impossibility of 
considering a mass assignment to this set ({12]). The mass of the set Uis then 
assigned to the empty set ø. Thus, the hybrid combination is defined ([13]) by 
Equations (5)-(8): 


m(A)=¢(A)*| H,(A)+H,(A)+H,(A)], VAeD? (5) 
H,(A)= > m,(X,)*m, (X,)*---*m, (X, ) (6) 
X; A, 
H,(A)= £ m (X,)*m,(X,)*---*m, (X,) (7) 
X,(i=1--,n)egd 


[i u(x; v-a (0 u(x; je#}n(a=0)| 
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H,(A)= m,(X,)*m, (X,)*- 
X;(i=1-,n)eD? 
X{UX_U UX, =A 
X1OX20 OX n=% 


«m,(X,) (8) 


The function ¢(A) of Equation (5) is a binary function equal to zero for empty 
or impossible sets A and is worth the unit for the others. In Equation (7), the set 
u(X;) represents the union of all the objects of the set X, Thus, this equation in- 
dicates that if the union of objects is also a constraint, then the mass is either as- 
signed to the union of all the singletons that form the objects, or to total ignor- 
ance. Equation (8) indicates that the mass is attributed to the union of all the 
objects of the sets X, if the intersection is a constraint. 

Several rules exist to achieve a hybrid combination. The 5th version of the fam- 
ily of rules of combination with Proportional Redistribution of Conflict (PCRS5) 
is presented in the following, because it gives better results ([14]). 

In a simplifying process, Djiknavorian ([12]) proposes a dynamic procedure 
which can be put in the form of the following algorithm (Dj): 

Beginning 
sı =m, 
For i = 1 to n—1 do 


Siy (C) = > Si (A)*m,,, (B) 


A,BeD 
AQNB=C 





DeD? S; (c) FM D) Mis 
DAC=0 

End 

m(C) =S, (c), YC eD’ 

End 


where m, (i = 1 to n) represent n mass functions corresponding to n different 
information sources and m the combined mass function obtained by the PCR5 
rule. 

The most commonly used decision rules for DSmT are based on combined 
mass functions, credibility functions, plausibility functions or pignistic probabil- 
ity functions. 

The functions of credibility (Bel), of plausibility (P/s) and of pignistic proba- 
bility (BetP) are defined from D? in [0,1] and are given respectively by Equa- 
tions (9)-(12): 


Bel(A)= > m(B), VAeD® (9) 

BeD? ,BcA 
Pls(A)= m(B), VAeD* (10) 

BeD? ,BOA+Ø 
BetP(A)= > BetP(6,), VA €D’ (11) 

QEA 
where 
m(A) . 

BetP(8.)= X —~, V@eO i=lp-,n (12) 


> 
AeD? b;eA |A| 
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|A] is the cardinal of the set A 

Moreover, the interval | Bel (A), Pls (A) | , known as the confidence interval, 
which quantifies the ignorance of the source on the hypothesis A, can be used 
for a decision. 

Following this presentation of the DSmT, the material used in this study is 


presented in the next section. 


3. Materials 


The tools used are software and data. 

With regard to the software, it was first used ENVI 4.7 to preprocessing ASTER 
images, then MATLAB to develop a model based on the use of the spectral in- 
dices NDVI, MNDWI and NDBa, and the DSmT to characterize states of the ve- 
getation, aquatic and mineral surfaces. The data are of two types: field data and 
remote sensing data. 

Field data consist of geographical coordinates of fixed points and outcrops. 
They are collected in the Toumodi-Yamoussoukro-Tiébissou area in the V Baoulé 
in Côte d'Ivoire. The remote sensing data used are derived from the ASTER sen- 
sor and are rectified satellite images of the scene AST_L1A_00301102004105832, 
covering this area. This sensor has 14 bands with a broad spectral region covering 
the visible and near infrared (VNIR-Visible and Near Infrared), the medium infra- 
red (SWIR-Short-Wave Infrared: Tape 4, Band 5, Band 6, Band 7, Band 8 and Band 
9) and Thermal Infrared (TIR-Thermal Infrared: Band 10, Band 11, Band 12, Band 
13 and Band 14). 

The spatial resolution associated with the said images is 15 m in the visible 
and the near infrared, 30 m in the medium infrared and 90 m in the thermal 


infrared. 


4. Methods 


The approach used consisted first of a preprocessing on the ASTER satellite im- 
ages under ENVI, and then it was developed a classification model based on the 
calculation of spectral indices (NDVI, MNDWI and NDBal) and the use of 
DSmT. Concretely, it was a question of modeling the discernment and reasoning 
frameworks, the mass functions as well as the functions of measuring the evi- 
dence, and defining the decision criteria. In addition, algorithms and program- 
ming codes in language were realized under Matlab software and the classified 


image was generated and evaluated. 


4.1. Preprocessing 


In order to benefit from the totality and the quality of the spatial resolutions and 
the spectral resolutions, the said ASTER satellite images have been subject to geo- 
referencing, geometric correction and resampling to create a compatible database, 
from the 14 bands. 


First, georeferencing was performed for each band using the k-nearest neigh- 
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bors method; then the geometric correction was made from 100 bitter points, 
chosen covering uniformly the ASTER scene of interest, with the bilinear me- 
thod; finally, the sampling, at a step of 15 m with the bilinear method, is carried 
out for the SWIR (bands 4, 5, 6, 7, 8 and 9) and TIR (bands 10, 11, 12, 13 and 14) 
bands. 

Georeferencing and geometric correction make it possible to make these satel- 
lite images superimposable on others georeferenced supports in the same coor- 


dinate system. 


4.2. Development of the Model 


4.2.1. Modeling of Discernment and Reasoning Frameworks 
Any portion of the Earth’s surface can be a combination of three main entities: a 
vegetation surface, an aquatic surface and a mineral surface. 

In this study, a vegetation area is an area of natural and/or cultural plants; 
an aquatic surface is a zone of natural and/or artificial watercourses and/or water 
bodies; a mineral surface is an area covered by soil, rock outcrops and/or built-up. 
Thus, depending on the size of the area of any surface observed, it is possible to 
physically observe a vegetation surface, an aquatic surface, a mineral surface, a ve- 
getation and aquatic surface, a vegetation and mineral surface, a mineral and aq- 
uatic surface or a vegetation, aquatic and mineral surface. The smaller the sur- 
face portion, the less it will contain different entities. So, an area of 15 m x 15 m 
could discriminate, at best on a satellite image according to the spectral charac- 
teristics, vegetation surfaces, aquatic surfaces and mineral surfaces. 


Therefore, the adopted discernment framework © is (Equation (13)): 
©={E,V,M} (13) 


where 
E: aquatic surface 
V: vegetation surface 
M: mineral surface 
By setting 0, = Æ 0, = V, 0, = Mand AUB={A,B} VA,B € D°, the follow- 
ing reasoning framework ([1]) is obtained: 
D° ={¢,E,V, M, {E,V}, {V,M}, {E,M}, ©, EAV,VAM,EOM, 
EQV,M},MO{EV},VO{E,M},{V,EOM},{M, EOV}, 
{EV AM}, [ENM,EAV,V AM},ENAV AM} 


4.2.2. Modeling Information Sources 

The sources of information considered in this study are the images produced 
by the neocanals obtained from the calculation of the spectral indices NDVI, 
MNDWIand NDBal. These sources have been chosen because of their simplici- 
ty of calculation, their standardized character and their reputation for less sensi- 
tivity (compared to reflectances) with respect to external factors such as the geo- 


metry of the illumination or the atmospheric effects. The calculation of the spec- 





DOI: 10.4236/ojapps.2017.710038 


526 Open Journal of Applied Sciences 


J.-C. Okaingni et al. 


tral indices NDVI, MNDWI and NDBal is obtained respectively by Equations 
(14)-(16): 











NDE p(bande 3) — p (bande 2) a4) 
p(bande 3) + p(bande 2) 
MNDWI = p(bande 1) — p(bande 7) (15) 
p(bande 1) + p(bande 7) 
NDBal = p(bande 14) - p(bande 7) (16) 
p (bande 14) + p(bande 7) 


p (bande 1): reflectance in the green (visible); 

p(bande 2) : reflectance in the red (visible); 

fe) (bande 3) : reflectance in the near infrared; 

p(bande 7) : reflectance in the mean infrared; 

p(bande 14) : reflectance in thermal infrared. 

The values of these indices are theoretically between —1 and +1. For NDVI, 
negative values correspond to surfaces such as snow, water or clouds, for which 
the reflectance in the red is greater than that of the near infrared. For bare soils, 
with reflectance values roughly the same in red and near infrared, the NDVI has 
values close to 0. Vegetation formations have positive NDVI values, generally 
between 0.1 and 0.7. The highest values correspond to the densest vegetation 
cover. The MNDW/I clearly distinguishes water surfaces (clear colors: positive val- 
ues) from non-water surfaces (vegetation, castings, bare soil) (dark colors: nega- 
tive values) giving more details than the NDWI ([15]). As for NDBal, it discri- 
minates the mineral surfaces (light colors: positive values) of the other compo- 
nents (dark colors: negative values) ([16]). 

For a successful work, a detection of segmentation thresholds was carried out 
by learning for each source by reviewing those obtained by the said authors. 
Thus, the thresholds used and retained are recorded in the Table 1. 


4.2.3. Modeling and Estimation of Masse Functions 
The mass functions of the sources are defined on the reasoning framework D? 


D? ={9,E,V,M,{E,V},{V,M}, {E,M}, ©, EQV,VOM,EOM, 
EQV,M\,MO{E,V},VO{E,M},{V,EOM},{M, EOV}, 
{E,V OM}, {EAM,EOV,V OAM}, EAV AM} 


Considering the normal distribution of variable x and parameters uA and dA, in 
Equation (17): 


Table 1. Segmentation thresholds considered for NDVI, MNDWIand NDBal. 





NDVI= X MNDWI= Y NDBal= Z 





Threshold X<-09 -09<X<0.1 01<X Y<S09 O9<Y 2Z<-01 -O01<Z 


Entity E M V {mV} E {BV M 
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(17) 


ea) 


1 
N (x, M4204) = > l TE 
A A 


with u, et a, respectively the mean and the standard deviation of the data x be- 
longing to A, the mass functions of the sources are then defined by Equations 
(18)-(24). 
- NDVI Function Mass 
With NDVI (x) : value of the pixel x of the NDVI image, we have: 
if NDVI (x) <-0.9 then: 
NDVI, (E)=0,V2nN (NDVI (x, Up,2,)) 
NDVI, (©) =1- NDVI, (E) (18) 
NDVI,(A)=0 WVAeD* -{E,} 
if —0.9 < NDVI (x) <0.1 then: 

NDVI, (M ) = oy V2nN (NDVI (x), fy +O ) 
NDVI,(©)=1-NDVI,(M) (19) 
NDVI,(A)=0 VAeD*-{M,0} 
if NDVI (x) >0.1 then: 

NDVI, (V)=o,V2nN (NDVI (x), 4,0, ) 

NDVI,(@)=1- NDVI, (V) (20) 

NDVI,(A)=0 YAeD’-{V,0} 

- MNDWI Function Mass 


With MNDWI ( x) : value of the pixel x of the MNDWTimage, we have: 
if MNDWI(x)<0.9 then: 


MNDWI, ({M,,V}) = Cu v} VEN (MNDWI (x), tiu v0 u v ) 
MNDWI,(©)=1- MNDWI,({M,V}) (21) 
MNDWI,(A)=0 WAeD°-{{M,V},0} 


if MNDWI (x) >0.9 then: 
MNDWI,(E) = 0,V2nN(MNDWI (x), u,,0,) 
MNDWI, (©) =1-MNDWI,(E) (22) 
MNDWI,(A)=0 WVAeD* -{E,@} 
- NDBal Function Mass 
With NDBal (x) : value of the pixel x of the NDBal image, we have: 
if NDBal (x) <-0.1 then: 
NDBal,, ({E,V}) = py) N2"N (NDBal (x), Hey)>F ev) ) 
NDBal (©) =1— NDBal, ({E,V}) (23) 
NDBal,(A)=0 VAeD*-{{E,V},0} 
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if NDBal (x) >-0.1 then: 


NDBal, (M )= cy V2TN (NDBal (x), fy.) 
NDBal (®)=1- NDBal,(M) (24) 
NDBal,(A)=0 WAeD*-{M,@} 


Focal elements of the sources S, = NDVI , S, =MNDWI and S, = NDBal are 

reported in Table 2. 

- Combined Mass Function 

The combined mass function is obtained by a hybrid combination method of 
DSmT based on the PCRS rule (algorithm (DJ), as a function of the twelve situ- 
ations generated by the thresholding conditions of S,(x), S,(x) and S,(x). 

It was used on the elements of D® and integrity constraint that generates four 
models that can reflect the reality on the ground: 

- model 1: integrity constraint: EOV OM =4¢; 

- model 2: integrity constraint! EAV AM =@ and cancellation of masses of 
EQ{V,M}, MO{E,V} and VO{E,M}; 

- model 3: integrity constraint! EAV AM =@ and cancellation of masses of 
EQ(V,M}, MOX{E,V}, VO{E,M}, {E,V}, {V,M} and {E,M}; 

- model 4: integrity constraint: EMV AM =@ and cancellation of masses of 
EQWV,M}, MOX{E,V}, VO{E,M}, {EV}, {V,M}, {EM}, 
EQV=¢, VAM=@¢ and EQOM=4¢; 

Thus, for each of the twelve situations, the combined mass function is gener- 
ated in the planes P,, P, and P, (P,P,P;), starting from the triplet of intersections 


formed by the focal elements of the sources S,, S, and S,, where: 


Table 2. Focal elements of the sources S,, S, and S, according to the thresholding condi- 
tionson S,(x), S,(x) et S,(x). 











Situation Conditions Focal elements 
al SWX WWR SWA S S, S; 
1 X, < -0.9 X,>0.9 X;2-0.1 {Ee} {E0} {M0} 
2 X, < -0.9 X,>09 X,<-0.1 {£0} {E0} {EV} 
3 X < -0.9 X% <09 X 2-0.1 {Be} {H{M, V},O} {M, ©} 
4 X; < -0.9 X% <0.9 X <-0.1 {Be} {{14V4,0} {BV} 
5 X >0.1 X>09 X%2-0.1 {V0} {EQ} {M0} 
6 X > 0.1 X%>09 X<-0.1 {VO} {£0} HE V},0} 
7 X, > 0.1 X% <09 X,>-0.1 {V0} {M 4,0} {M, ©} 
8 X > 0.1 X% <09 %<-01 {V0} HM V},0} HA v0} 
9 -0.9<X,<01 X%>09 X% 2-01 {MO} {EQ} {Me} 
10 -0.9<X,<01 %>09 X;,-0.1 {M0} {£0} {{B n0} 
11 -0.9<X,<0.1 X,<0.9 X,2-0.1 {M,0} {{M,V},0} {M,@} 
12 -0.9<X,<0.1 X,<09 X;<-01 {MO} {{M,V},} {{E,V},0} 
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P, ={E,V,M,@} = {1;2;3;4} 
P, ={E,{M,V},@} ={1;2;3} 
P,={M.{E,V},0} = {12:3} 


In the planes P,P,P;, the first component of the intersection belongs to the 
plane P, and is indicated by its position in the same plane. The second and third 
components obey the same principle respectively in the planes P, and P}. 


4.2.4. Decision and Evaluation 

Once all the combined mass functions of the single and multiple assumptions of 
a pixel x are determined, the decision rule used is the combined mass function 
maximum given the presence of intersections of elements in the model of DSmT 
({1)). 

The evaluation for this study is based first on statistical analysis and then on a 
visual conformance analysis. The statistical analysis is carried out by a confusion 
matrix Mok) based on the ground truth classes and those of the combined im- 
age by model k (k= 1, 2, 3, 4). True field classes are placed in columns while 
those in the combined image are in rows. The total number of pixels per class for 
the ground truth is distributed to the classes of the combined image. The perfor- 


mance indices are then calculated from Equations (25)-(28) ([2] [17]): 


oh fas 
cor = Ma) i (25) 
N, 
3 
X GCR, 
GCR roy E m (26) 
M,y(ii) . 
uea T (27) 
Mta 
gcr = y Me) (28) 
jji N; 


with # class number of truth field, 

$ class number of the combined image by model, 

N; number of pixels of class of the ground truth, 

Mji): number of pixels of the class / of the ground truth having been as- 
signed after classification to the class j of the combined image by model, 

GCR; rate of well classified pixels of class of ground truth, 

GCR moy: average rate of well classified pixels of truth ground, 

ECR;; rate of misclassified pixels of class j of ground truth in class j of the 
combined image by model, 

ECR; rate of misclassified pixels of class j of the ground truth. 

Finally, a synthetic confusion matrix Mc, is produced by combining the four 
models. 


The ground truth was realized from a field reconnaissance to highlight the enti- 
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ties water “EZ”, vegetation “ V” and mineral “M”. Once these entities are found, from 
their geographical coordinates, they are identified on the images and then sam- 
pled using the pixels. 

The visual conformance analysis, on the other hand, consisted of verifying on 
the ground the correspondences of the various composite classes provided by the 
classification. Portions of the image have been chosen and their geographical 
coordinates determined, with which a field check is carried out. 


The methodological approach used is summarized in Figure 1. 


The different results obtained during this process are presented in the follow- 


Data : ASTER Images 
scene AST_L1A_00301102004105832 


Preprocessing : 
- Georeferencing 


ing section. 


- Resampling 
- Geometric correction 
- Data base creation 


Data base: 
Bands: from 1 to 14 
Spectral indices computation 
NDVI segmented (E, V , M) MNDWI segmented (E, {V,M}) 


DSmT Modeling 
NDVI mass function MNDWI mass function 
defined on D? defined on D® 
DSmT by Hybrid Combination 
with Rule PCRS 
Combined mass function 
defined on D? 
Decision making by maximum of 
combined mass function 


Combined Image 
defined on D® 
















Segmentation 











NDBal segmented ({E,V}, M) 


NDBal mass function 
defined on D® 
Integrity constraints 



















Figure 1. The flow chart of the methodology used. 
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5. Results and Discussions 


The results presented derive from the classifications made with the hybrid DSmT 
based on the PCRS rule for the above four models generated by an integrity con- 


straint and simplifications. 


5.1. Classification Models 


Under the constraint of exclusion integrity EAV AM =¢ for the model 1 and 
neglecting the masses of EN{V,M}, MA{E,V}, VA{E,M}, {EV}, {VM}, 
{E,M}, EAV, MAV and EAM in accordance with models 2, 3 and 4, 
Figures 2-5 respectively show the results of the classification of models 1, 2, 3 
and 4 whose number of classes produced according to said models is 11, 8, 6 and 
3. From a qualitative point of view, according to the model, we observe three 
simple classes for certain knowledge (£, V, M) and/or eight composite classes. 
The latter are formed by two classes of union characterizing partial ignorances 
({E,V}, {M,V}) and six intersection classes including three simple intersec- 
tion classes (EOV, MOV, ECM ) and three classes of composite intersec- 
tion (EN{M,V}, MO{E,V}, VO{E,M}), which represent paradoxes. The 
distribution of pixels for quantitative analysis in these different classes of images 
classified according to models 1, 2, 3 and 4 is recorded in Table 3. In modeling, 
model 1 was retained with 84.98% for single classes versus 15.02% for composite 
classes. 


5.2. Evaluation 


The evaluation was done successively from models 4, 3, 2 and 1. The statistical 
analysis is presented in Table 4. From model 4, well classified pixel rates for E 
(96.02%), V (91.51%) and M (99.60%) are obtained, giving an average rate of 
95.74%. The pixels E of the ground truth are coincident with the classification 
pixels V and M respectively 0.29% and 3.69%, whereas the pixels V and M of the 
ground truth are mutually coincident with 8.41 % and 0.40%. For models 1, 2 
and 3, in addition to the statistical analysis, which gives a minimal percentage of 





Figure 2. Image classified according to model 1. 
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Figure 4. Image classified according to model 3. 
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Figure 5. Image classified according to model 4. 
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Table 3. Distribution of pixels by class according to models 1, 2, 3 and 4. 
































Model 1 Model 2 Model 3 Model 4 
pisces % pixels saath % pixels nares % pixels Presid % pixels 
E 129,600 0.75 226,870 1.32 289,447 1.68 369,509 2.15 
V 4,767,345 27.68 4,967,525 28.85 5,250,117 30.49 5,288,132 30.71 
M 9,737,965 56.55 10,562,346 61.34 11,439,968 66.43 11,562,359 67.14 
{E V} 208,329 1.21 566,710 3.29 0 0 0 0 
{M, V} 553,069 3.21 771,473 4.48 0 0 0 0 
Combined 
and {EM} 0 0 0 0 0 0 0 0 
classified EnV 29,967 0.17 29,967 0.17 29,967 0.17 0 0 
image 
LIMAV 1714 0.01 1714 0.01 75,627 0.44 0 0 
EnM 68,047 0.40 93,395 0.54 134,874 0.78 0 0 
EN{M, V} 218,685 1.27 0 0 0 0 0 0 
MANE, V} 1,092,685 6.35 0 0 0 0 0 0 
WA{ MF} 412,594 2.40 0 0 0 0 0 0 
TOTAL 17,220,000 100 17,220,000 100 17,220,000 100 17,220,000 100 
Table 4. Synthesis confusion matrix Mcp combining the four models. 
Truth on the ground 
Model 1 Model 2 Model 3 Model 4 
E Vv M E vV M E vV M E vV M 
E 94.24 0.00 0.00 94.30 0.00 0.00 94.44 0.00 0.00 96.02 0.00 0.00 
Vv 0.00 90.67 0.00 0.00 91.19 0.00 0.00 91.27 0.00 0.29 91.59 0.40 
M 0.00 0.11 95.11 0.00 0.17 95.63 0.23 6.41 99.58 3.69 8.41 99.60 
{E,V} 0.00 3.14 0.00 0.00 3.14 1.25 0.00 0.00 0.00 0.00 0.00 0.00 
{M, V} 0.00 4.14 0.00 0.00 4.49 0.02 0.00 0.00 0.00 0.00 0.00 0.00 
Combined 
and {E,M} 0.62 0.00 0.04 0.70 0.00 0.16 0.00 0.00 0.00 0.00 0.00 0.00 
classified gay 0.00 091 000 0.04 0.91 0.00 000 0.99 000 000 000 0.00 
SAES LMoV 0.00 0.00 1.97 0.00 0.09 1.97 0.00 1.32 0.42 0.00 0.00 0.00 
EQ M 4.79 0.00 0.96 4.95 0.00 0.96 5.33 0.00 0.00 0.00 0.00 0.00 
En{M, V} 0.35 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
MANE, V} 0.00 0.00 1.92 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
VA{ M, E 0.00 1.02 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 000 
ECR; 5.76 9.33 4.89 5.70 8.81 4.37 5.56 8.73 0.42 3.98 8.41 0.40 





well classified pixels of 90.67%, the conformity analysis carried out in the field 
gives 96.37%, 94.19% and 91.02%, an average rate of 93.86%. Therefore, the ini- 


tial model 1 is chosen as a model adapted to this situation. 
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6. Conclusion 


In this article, in order to map the surface state of the Toumodi-Yamoussoukro- 
Tiébissou zone in the “V Baoule” of Côte d'Ivoire, the theory of plausible and 
paradoxical reasoning of Dézert-Smarandache was proposed with the aim of tak- 
ing into account the paradoxical nature of the intersections of the vegetation sur- 
face, the aquatic surface and the mineral surface from Aster satellite images and 
the NDVI, MNDWI and NDBal spectral indices. After a preprocessing on the 
ASTER satellite images, a model of surface state characterization based on the 
calculation of spectral indices (NDVI, MNDWIand NDBal) and the DSmT with 
the computation PCR5 rule, was developed. It gives an average rate of 93.34% 
for well classified pixels and a field compliance rate of 96.37%. However, the re- 
sults could be improved by a supervised approach to mass functions, in particu- 
lar through the extensive use of learning elements. This model could be used, 


with appropriate adjustments, for other mapping purposes. 
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